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Abstract 

The cross section for ep — > e bbX in photoproduction is measured with the HI detector 
at the ep-collider HERA. The decay channel bb — > eeX' is selected by identifying the semi- 
electronic decays of the 6-quarks. The total production cross section is measured in the 
kinematic range given by the photon virtuality Q 2 < lGeV 2 , the inelasticity 0.05 < y < 
0.65 and the pseudorapidity of the 6-quarks |r/(6)|, \rj(b)\ < 2. The differential production 
cross section is measured as a function of the average transverse momentum of the beauty 
quarks {Px{b)) down to the threshold. The results are compared to next-to-leading-order 
QCD predictions. 



Submitted to Eur. Phys. J. C 



F. D. Aaron 5,45 , C. Alexa 5 , V. Andreev 25 , S. Backovic 30 , A. Baghdasaryan 38 , 

S. Baghdasaryan 38 , E. Barrelet 29 , W. Bartel 11 , K. Begzsuren 35 , A. Belousov 25 , P. Belov 11 , 
J.C. Bizot 27 , V. Boudry 28 , 1. Bozovic-Jelisavcic 2 , J. Bracinik 3 , G. Brandt 11 , M. Brinkmann 11 , 
V. Brisson 27 , D. Britzger 11 , D. Bruncko 16 , A. Bunyatyan 13 ' 38 , A. Bylinkin 24 , L. Bystritskaya 24 , 
A.J. Campbell 11 , K.B. Cantun Avila 22 , F. Ceccopieri 4 , K. Cerny 32 , V. Cerny 16 , V. Chekelian 26 , 
J.G. Contreras 22 , J.A. Coughlan 6 , J. Cvach 31 , J.B. Dainton 18 , K. Daum 37 ' 42 , B. Delcourt 27 , 
J. Delvax 4 , E.A. De Wolf 4 , C. Diaconu 21 , M. Dobre 12 ' 47 ' 48 , V. Dodonov 13 , A. Dossanov 12 ' 26 , 
A. Dubak 30 , G. Eckerlin 11 , S. Egli 36 , A. Eliseev 25 , E. Elsen 11 , L. Favart 4 , A. Fedotov 24 , 
R. Feist 11 , J. Feltesse 10 , J. Ferencei 16 , D.-J. Fischer 11 , M. Fleischer 11 , A. Fomenko 25 , 

E. Gabathuler 18 , J. Gayler 11 , S. Ghazaryan 11 , A. Glazov 11 , L. Goerlich 7 , N. Gogitidze 25 , 
M. Gouzevitch 11 - 43 , C. Grab 40 , A. Grebenyuk 11 , T. Greenshaw 18 , G. Grindhammer 26 , 

S. Habib 11 , D. Haidt 11 , R.C.W. Henderson 17 , E. Hennekemper 15 , H. Henschel 39 , M. Herbst 15 , 

G. Herrera 23 , M. Hildebrandt 36 , K.H. Hiller 39 , D. Hoffmann 21 , R. Horisberger 36 , T. Hreus 4 , 

F. Huber 14 , M. Jacquet 27 , X. Janssen 4 , L. Jonsson 20 , A.W. Jung 15 ' 51 , H. Jung 11 ' 4 , 

M. Kapichine 9 , 1.R. Kenyon 3 , C. Kiesling 26 , M. Klein 18 , C. Kleinwort 11 , R. Kogler 12 , 
R Kostka 39 , M. Kramer 11 , J. Kretzschmar 18 , K. Kriiger 15 , M.RJ. Landon 19 , W. Lange 39 , 

G. Lastovicka-Medin 30 , R Laycock 18 , A. Lebedev 25 , V. Lendermann 15 , S. Levonian 11 , 

K. Lipka 11,47 , B. List 11 , J. List 11 , B. Lobodzinski 11 , R. Lopez-Fernandez 23 , V. Lubimov 24 , 

E. Malinovski 25 , H.-U. Martyn 1 , S.J. Maxfield 18 , A. Mehta 18 , A.B. Meyer 11 , H. Meyer 37 , 

J. Meyer 11 , S. Mikocki 7 , 1. Milcewicz-Mika 7 , F. Moreau 28 , A. Morozov 9 , J.V. Morris 6 , 

K. Miiller 41 , Th. Naumann 39 , RR. Newman 3 , C. Niebuhr 11 , D. Nikitin 9 , G. Nowak 7 , 

K. Nowak 12 , B. Olivier 26 , J.E. Olsson 11 , D. Ozerov 11 , P. Pahl 11 , V. Palichik 9 , M. Pandurovic 2 , 

C. Pascaud 27 , G.D. Patel 18 , E. Perez 10 ' 44 , A. Petrukhin 11 , 1. Picuric 30 , H. Pirumov 14 , D. Pitzl 11 , 
R. Placakyte 11 , B. Pokorny 32 , R. Polifka 32 ' 49 , B. Povh 13 , V. Radescu 11 , N. Raicevic 30 , 

T. Ravdandorj 35 , P. Reimer 31 , E. Rizvi 19 , P. Robmann 41 , R. Roosen 4 , A. Rostovtsev 24 , 
M. Rotaru 5 , J.E. Ruiz Tabasco 22 , S. Rusakov 25 , D. Salek 32 , D.P.C. Sankey 6 , M. Sauter 14 , 

E. Sauvan 21 ' 50 , S. Schmitt 11 , L. Schoeffel 10 , A. Schoning 14 , H.-C. Schultz-Coulon 15 , 

F. Sefkow 11 , L.N. Shtarkov 25 , S. Shushkevich 11 , T. Sloan 17 , Y. Soloviev 11 ' 25 , P. Sopicki 7 , 

D. South 11 , V. Spaskov 9 , A. Specka 28 , Z. Staykova 4 , M. Steder 11 , B. Stella 33 , G. Stoicea 5 , 
U. Straumann 41 , T. Sykora 4 ' 32 , PD. Thompson 3 , TH. Tran 27 , D. Traynor 19 , P. Truol 41 , 

I. Tsakov 34 , B. Tseepeldorj 35 ' 46 , J. Tumau 7 , A. Valkarova 32 , C. Vallee 21 , P. Van Mechelen 4 , 
Y. Vazdik 25 , D. Wegener 8 , E. Wiinsch 11 , J. Zacek 32 , J. Zalesak 31 , Z. Zhang 27 , A. Zhokin 24 , 
R. Zlebcfk 32 , H. Zohrabyan 38 , and F. Zomer 27 

1 /. Physikalisches Institut der RWTH, Aachen, Germany 

2 Vinca Institute of Nuclear Sciences, University of Belgrade, 1 100 Belgrade, Serbia 

3 School of Physics and Astronomy, University of Birmingham, Birmingham, Uf^ 

4 Inter-University Institute for High Energies ULB-VUB, Brussels and Universiteit Antwerpen, 
Antwerpen, Belgium 

5 National Institute for Physics and Nuclear Engineering (NIPNE) , Bucharest, Romania k 

6 STFC, Rutherford Appleton Laboratory, Didcot, Oxfordshire, UK b 

7 Institute for Nuclear Physics, Cracow, Poland d 

8 Institut fur Physik, TU Dortmund, Dortmund, Germany 11 

9 Joint Institute for Nuclear Research, Dubna, Russia 

10 CEA, DSM/Irfu, CE-Saclay, Gif-sur-Yvette, France 

11 DESY, Hamburg, Germany 



1 



12 InstitutfUr Experimentalphysik, Universitat Hamburg, Hamburg, Germany 11 

13 Max-Planck-Institut ftir Kernphysik, Heidelberg, Germany 

14 Physikalisches Institut, Universitat Heidelberg, Heidelberg, Germany 11 

15 Kirchhoff- Institut fur Physik, Universitat Heidelberg, Heidelberg, Germany a 

16 Institute of Experimental Physics, Slovak Academy of Sciences, Kosice, Slovak Republic e 

17 Department of Physics, University of Lancaster, Lancaster, UK b 

18 Department of Physics, University of Liverpool, Liverpool, UK 1 * 

19 School of Physics and Astronomy, Queen Mary, University of London, London, UK b 

20 Physics Department, University of Lund, Lund, Sweden? 

21 CPPM, Aix-Marseille Univ, CNRS/IN2P3, 13288 Marseille, France 

22 Departamento de Fisica Aplicada, CINVESTAV, Merida, Yucatan, Mexico 1 

23 Departamento de Fisica, CINVESTAV IPN, Mexico City, Mexico 1 

24 Institute for Theoretical and Experimental Physics, Moscow, Russia^ 

25 Lebedev Physical Institute, Moscow, Russia 

26 Max-Planck-Institut fur Physik, Munchen, Germany 

27 LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 

28 LLR, Ecole Polytechnique, CNRS/IN2P3, Palaiseau, France 

29 LPNHE, Universite Pierre et Marie Curie Paris 6, Universite Denis Diderot Paris 7, 
CNRS/IN2P3, Paris, France 

30 Faculty of Science, University of Montenegro, Podgorica, Montenegro 1 

31 Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech Republic 9 

32 Faculty of Mathematics and Physics, Charles University, Praha, Czech Republic 9 

33 Dipartimento di Fisica Universita di Roma Tre and INFN Roma 3, Roma, Italy 

34 Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria 

35 Institute of Physics and Technology of the Mongolian Academy of Sciences, Ulaanbaatar, 
Mongolia 

36 Paul Scherrer Institut, Villigen, Switzerland 

37 Fachbereich C, Universitat Wuppertal, Wuppertal, Germany 

38 Yerevan Physics Institute, Yerevan, Armenia 

39 DESY, Zeuthen, Germany 

40 InstitutfUr Teilchenphysik, ETH, Zurich, Switzerland 11 

41 Physik-Institut der Universitat Zurich, Zurich, Switzerland 11 

42 Also at Rechenzentrum, Universitat Wuppertal, Wuppertal, Germany 

43 Also at IPNL, Universite Claude Bernard Lyon 1, CNRS/IN2P3, Villeurbanne, France 

44 Also at CERN, Geneva, Switzerland 

45 Also at Faculty of Physics, University of Bucharest, Bucharest, Romania 

46 Also at Ulaanbaatar University, Ulaanbaatar, Mongolia 

47 Supported by the Initiative and Networking Fund of the Helmholtz Association (HGF) under 
the contract VH-NG-401. 

48 Absent on leave from NIPNE-HH, Bucharest, Romania 

49 Also at Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 

50 Also at LAPP, Universite de Savoie, CNRS/IN2P3, Annecy-le-Vieux, France 

51 Now at Fermi National Accelerator Laboratory, Batavia, Illinois 60510, USA 

a Supported by the Bundesministeriumfur Bildung und Forschung, FRG, under contract 
numbers 05H09GUF, 05H09VHC, 05H09VHF, 05H16PEA 



2 



b Supported by the UK Science and Technology Facilities Council, and formerly by the UK 
Particle Physics and Astronomy Research Council 

c Supported by FNRS-FWO-Vlaanderen, IISN-IIKW and IWT and by Interuniversity Attraction 
Poles Programme, Belgian Science Policy 

d Partially Supported by Polish Ministry of Science and Higher Education, grant 

DPN/N168/DESY/2009 

e Supported by VEGA SR grant no. 2/7062/27 

* Supported by the Swedish Natural Science Research Council 

9 Supported by the Ministry of Education of the Czech Republic under the projects LC527, 

INGO-LA09042 and MSM002 1620859 

h Supported by the Swiss National Science Foundation 

1 Supported by CONACYT, Mexico, grant 48778-F 

i Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2 and Rosatom 

k Supported by the Romanian National Authority for Scientific Research under the contract PN 

09370101 

1 Partially Supported by Ministry of Science of Montenegro, no. 05-1/3-3352 



3 



1 Introduction 



In ep collisions at HERA beauty quarks are mainly produced as bb pairs via the fusion of a 
quasi-real photon emitted by the incoming electron (or positron) and a gluon of the proton 
as depicted in figure [TJl. This process is referred to as direct or pointlike and can be calculated 
using perturbative quantum chromodynamics (QCD) due to the large scale provided by the mass 
of the heavy 6-quark and the correspondingly small coupling a s . Resolved processes where 
the photon fluctuates into a hadronic state before undergoing a hard collision, as indicated in 
figure [Tt>, are expected to be largely suppressed compared to the direct production process, 
because of the large 6-quark mass. Due to the dominance of the direct process over the resolved 
process, the production of 6-quarks in ep collisions at HERA is an excellent testing ground for 
QCD predictions. 



a) b) 




Figure 1: Generic leading order diagrams for bb production in ep collisions. The diagram a) is 
referred to as direct or pointlike, the diagram b) is referred to as resolved or hadronlike. 

Theory uncertainties in the prediction of the cross section, which are mainly related to the 
renormalisation and factorisation scales, are expected to be smaller for beauty production than 
for charm production. The study of beauty photoproduction near threshold is of particular 
theoretical interest as the only hard scale in this process is provided by the 6-quark mass, and 
other scales like the photon virtuality (Q 2 fa GeV 2 in photoproduction) or the transverse 
momentum of the 6-quark can be neglected. 

At HERA the beauty cross section in photoproduction ep — > e bbX has been measured by 
the HI [Q]-[5l and ZEUS Il6l-fl2l collaborations and compared to calculations lfT3til"5l at next-to- 
leading order (NLO) QCD, performed in the fixed flavour number scheme in which the beauty 
quark is treated as massive. In general the predictions using the factorisation and renormali- 
sation scale fiR = /if — \/ m \ + Pr(b) 2 do not agree well with the data. In particular at low 
values of the transverse momentum of the beauty quarks Pr(b) ~ GeV, i.e. in the phase space 
region where the only hard scale involved is the fe-quark mass mi, the measurements show a ten- 
dency to lie above the prediction. The choice of a lower scale, fi R = fi F = 1/2 ^/m 2 + P T (b) 2 , 
leads to a better agreement of the prediction with the data f[T6l . 

In the present analysis a measurement of the differential beauty cross section at HERA in 
photoproduction as function of the quadratically averaged transverse momentum of the pro- 
duced beauty quarks, da/d(PT(b)), is made down to the ^-production threshold, using a novel 
technique based on low momentum electron identification. 
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Most of the previous beauty measurements at HERA in photoproduction and deep-inelastic 
scattering (DIS, Q 2 > 1 GeV 2 ) identified jets of 6-quarks using single leptons tags [Q3Sm[Z]|9]- 
[TTl[T7] - |20l or displaced vertices fl2ll6l l2TI - l24l . Jet-based 6-tagging algorithms are in general very 
efficient at high transverse momenta of 6-quarks, Px(b) > 6 — 7 GeV, but degrade significantly 
for lower values due to the absence of the boost and the short decay length. In addition 6-quarks 
almost at rest lead to isotropic decay topologies of the final state where jet finders usually 
fail. A second class of analyses used double tags to select bb pairs either by reconstructing 
two muons Jl or a muon and a D* meson (3l[T0l, utilizing the semileptonic decay channel 
b — > fiX' and the decay channel b — > D*X', respectively. Lower values of the transverse 
momentum of the 6-quarks become accessible by the use of lepton tags without requiring jets, 
where the minimum Pr(b) value is determined by the minimum transverse momentum cut on 
the lepton. For muons this cut is typically at Pt((j) ~ 2 GeV, and therefore too high to measure 
efficiently the production cross sections of 6-quarks near threshold. 

In the present analysis the differential beauty cross section is measured using electron pairs, 
exploiting the double- semileptonic decay bb — > eeX', with online and offline Pr(e) thresholds 
for the electron identification of about 1 GeV. The events were recorded by identifying low 
momentum electrons already online using a dedicated trigger, which recorded data in the year 
2007 with a corresponding integrated luminosity of 48.1 pb _1 . This low cut on the transverse 
electron momentum, Pr{e), improves not only the total acceptance but also makes the low 
Pr(b) phase space experimentally accessible. 

2 Monte Carlo Simulations and QCD Calculations 

The Monte Carlo generators PYTHIA [j25] and CASCADE [|26l are used to determine the signal 
efficiency and the detector acceptance for the process ep — > e bbX — > e eeX', and to simulate 
the production of charm quarks. Differences in the predictions are taken into account as system- 
atic uncertainty, see section [5~4l For the production of J/ip mesons only CASCADE and for the 
production of light quarks in photoproduction only PYTHIA is used. Deep-inelastic scattering 
is simulated using the Monte Carlo generator RAPGAP [|27l . 

In PYTHIA leading order matrix elements are implemented taking into account the mass 
of the heavy quarks. The CTEQ6L [28J set of proton parton density functions is used. The 
parton shower evolution in PYTHIA is based on the DGLAP equations [|29l . In addition to the 
direct process, the resolved photon component is calculated by using the photon parton density 
function SAS 2D QUI . 

For the CASCADE simulation the direct — > bb and 7*p — > cc processes are imple- 
mented using off-shell matrix elements, which are convoluted with fc r -unintegrate<£] proton 
parton density functions. The AO [|3T1 set of parton density functions is used. The parton evo- 
lution in CASCADE is based on the CCFM evolution equation ll32l for the initial state parton 
shower. 

In PYTHIA, CASCADE and RAPGAP higher order QCD corrections are included by sim- 
ulating parton showers in the initial and final state. These Monte Carlo generators use the Lund 

1 kx denotes the transverse momentum of the parton. 
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String Model 11331 for simulating the hadronisation of light quarks. For the hadronisation of 
heavy quarks the Bowler fragmentation model [1341 is employed with parameters as used in a 
previous analysis ETI . 

In order to correct for detector effects and to estimate the systematic uncertainties associated 
with the measurement, the generated events are passed through a detailed simulation of the 
detector response based on the GEANT program [35 1 and through the same reconstruction and 
analysis software as is used for the data. 

Theory cross sections are calculated in NLO QCD in the fixed flavour number scheme using 
the program FMNR [[T3Tfl5l in order to compare with the data. These calculations are expected 
to give reliable results in the kinematic region considered here, where the transverse momentum 
of the heavy quark is of the same order of magnitude as its mass. The calculations are performed 
as a function of the quadratically averaged transverse momentum of the produced beauty pair 

(Mb)) = . (1) 

The prediction of FMNR is evaluated for the direct and resolved photon processes. For the pro- 
ton the CTEQ6M ll28l set and for the photon the GRV-HO 11361 set of parton density functions 
are used. In this analysis, the renormalisation and factorisation scales are chosen to be equal, 
— fJ-F — A*o> with /i = 1/2 y/m\ + {Prib)) 2 and rrib = 4.75 GeV. The value used for the 
QCD scale Aq CD corresponds to the value of the strong coupling constant a s (M z ) =0.118. The 
theoretical uncertainty of the prediction is evaluated by varying the scales fiR and fip simultane- 
ously in the window /xo/2 < [Ir,f < 2/io and the beauty mass in the range 4.5 < m b < 5.0 GeV. 
By recalculating the cross section with different parton density functions the theoretical uncer- 
tainty due to the choice of the photon and proton parton density functions is found to be much 
smaller than the theoretical uncertainties and thus is neglected. 



3 HI Detector 

A detailed description of the HI detector can be found elsewhere [137113811 . In the following, only 
detector components relevant to this analysis are briefly discussed. The origin of the HI coordi- 
nate system is the nominal ep interaction point, with the direction of the proton beam defining 
the positive z axis (forward direction). Transverse momenta are measured in the x-y plane. 
Polar (#) and azimuthal (0) angles are measured with respect to this reference system. The 
pseudorapidity is defined to be rj = — In tan($/2). 

In the central region (15° < i3< 165°) the interaction point is surrounded by the central 
tracking detector (CTD). The CTD comprises two large cylindrical jet chambers (CJC1 and 
CJC2) and the silicon vertex detector [|39l . The CJCs are separated by a drift chamber which 
improves the z coordinate reconstruction. The CTD detectors are arranged concentrically 
around the interaction region in a solenoidal magnetic field of 1.16 T. The trajectories of 
the charged particles are measured with a transverse momentum resolution of <j(pt)/pt ~ 
0.2% pt /GeV © 1.5%. In addition the CJCs provide a measurement of the specific ionisation 
energy loss dE/dx of charged particles with a relative resolution of 6.5% for long tracks. A set 
of five cylindrical multiwire proportional chambers [40] mainly used for first level triggering are 
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situated inside the inner CJC1 covering the polar angular region 11° < -d < 169°. The forward 
tracking detector and the backward proportional chamber measure tracks of charged particles at 
smaller (7° < ■d < 25°) and larger (155° < $ < 175°) polar angles than the central tracker, respec- 
tively. 

The liquid argon (LAr) sampling calorimeter [41 j surrounds the tracking chambers and has 
a polar angle coverage of 4° <$< 154°. It consists of an inner electromagnetic section with 
lead absorbers and an outer hadronic section with steel absorbers. The LAr calorimeter is di- 
vided into eight wheels along the beam axis. The electromagnetic and the hadronic sections are 
highly segmented in the transverse and the longitudinal directions. Energies of electromagnetic 
showers are measured with a precision of a(E)/E = 12%/v/E/GeV © 1% and energies of 
hadronic showers with a{E)/E = 50%/ y/E/ GeV © 2%, as determined in test beam experi- 
ments Il42ll43il . In the backward region (153° <$< 178°), particle energies are measured by a 
lead- scintillating fibre calorimeter (SpaCal) [38] . 

The luminosity is determined from the rate of the elastic QED Compton process ep — > e jp, 
with the electron and the photon detected in the SpaCal calorimeter, and the rate of DIS events 
measured in the SpaCal calorimeter [|44l . 

For data collection a four level trigger system is employed, of which the first two levels 
are implemented in hardware. The first level trigger (LI) is based on various sub-detector 
components, which are combined and refined at the second level (L2). The third level (L3) is a 
software based trigger using combined LI and L2 trigger information from various subdetector 
components. Fully reconstructed events are subject to an additional selection at the software 
filter farm (L4). 

The data used for this measurement were recorded by the Fast Track Trigger (FTT) H45 11461 
which, based on hit information provided by the CJCs, reconstructs tracks with subsequently 
refined granularity at the first two trigger levels, first in the x-y plane at LI and then in three 
dimensions at L2. Of special importance is the third trigger level integrated in the FTT fl46l . 
which identifies low energy electrons (E > 1 GeV) [|47ll48l by combining FTT tracks with 
energy depositions reconstructed in the LAr calorimeter by the Jet Trigger (JT) [49 ]. 



4 Experimental Method 

The data sample used for this analysis was recorded in the year 2007, when positrons at an 
energy of 27.6 GeV collided with protons at 920 GeV, and when all trigger levels of the FTT 
and the JT were in operation. The recorded data corresponds to a total integrated luminosity of 
48.1 pb _1 . 

In this analysis the measurement of bb photoproduction is based on the identification and 
selection of two electrons in the LAr calorimeter at low transverse momentum -Pr(e) > 1 GeV 
to tag the semi-electronic decays of the 6-quarks. In about 2% of all 66-decays two electrons 
originate from the same 6-quark from the decay chain b — > c e~u e — > s e~e + v e v e . In about 4% 
of all 66-decays the two electrons originate from decays of different fe-quarks, where they are 
either produced directly in the semi-electronic 6-decays or in the subsequent semi-electronic 
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Subtrigger 


# ele. cand. 


Pr, ftt [GeV] 


p T ' jt [%] 

"t.ftt 


Acp [rad] 


Atf [rad] 


£ [pb- 1 ] 


1ow-Pt 


> 2 


> 1.2 


> 30 


< 0.30 


< 0.25 


25.1 


medium-P T 


> 1 


> 1.5 


> 50 


< 0.15 


< 0.20 


13.1 


high-Pr 


> 1 


> 2.0 


> 60 


< 0.20 


< 0.20 


33.5 



Table 1: L3-online cuts used to trigger electron candidates. Explanations to the cuts are given 
in the text. The last column contains the prescale corrected integrated luminosity of each sub- 
trigger. The medium- P r subtrigger was commissioned at a later stage. 



c-decays. The electrons can be either of opposite charge (combinations bb — > e + e~X and 
cc — > e + e~X) or of same charge (combinations be — >■ e~e~X and be — > e + e + X). These charge 
relations hold only in the case of no BB mixing. In the following all possible combinations 
including BB mixing are considered in order to discriminate bb decays against semi-electronic 
decays of cc events. Electron pairs from J / ip decays are distinguished from those from 6-decays 
by reconstructing their invariant mass. Misidentified electrons originating mainly from the light 
quark background are constrained by varying the cuts on the electron identification described 
in section [51 

4.1 Online Electron Identification 

Events containing several tracks and one or two electron candidates compatible with the signa- 
ture of semi-electronic 6-decays are triggered, using the FTT on the trigger levels LI to L3. On 
the first trigger level more than five tracks with transverse momentum thresholds in the range 
0.1 — 1.8 GeV are required. These high multiplicities are verified at the second trigger level, 
exploiting the higher track resolution available at this level. On the third trigger level the track 
information as determined by FTT-L2 is combined with the energy depositions as measured in 
the LAr calorimeter by the Jet Trigger [146T - I481 to identify electrons. Electron candidates are 
required to fulfill a geometrical track-cluster matching condition using the distance variables 
A$ = I^ftt — $jt| and Aip = |v?ftt — </?jt|- hi addition the transverse momentum Pt,ftt 
as measured with the FTT-L2 has to be compatible with the associated transverse energy Pt,jt 
measured in the LAr calorimeter by the JT. A lower cut on the quantity E Tj jt/Pt,ftt is used 
to discriminate electrons against hadrons, which deposit significantly less energy in the non- 
compensating LAr calorimeter. 

For this analysis three subtriggers are used, which have identical LI and similar L2 trigger 
conditions, but different conditions on L3 as summarized in table [Q The subtrigger with the 
lowest transverse momentum threshold of Pt,ftt > 1.2 GeV requires events with at least two 
electron candidates. The other two subtriggers select events with a minimum of one electron 
candidate with Pt,ftt thresholds of 1.5 and 2.0 GeV. The three data sets recorded by these 
FTT-JT based subtriggers cover an overlapping kinematic phase space, but correspond to dif- 
ferent integrated luminosities due to different trigger prescale factors. The three data sets are 
combined using a weighting method [47] to account for correlated triggers with prescales. The 
individual prescale corrected luminosities are also given in table [T] 
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4.2 Offline Electron Identification 



Electrons in the polar angle range of 20° < $(e) < 140° and with a transverse momentum of 
Pr(e) > lGeV, with -Pr(e) and $(e) measured from the electron track, are identified using 
energy depositions in the LAr calorimeter and specific ionisation loss measured in the CJCs. 
Two estimators, -D calo and D dE / dx , are defined to discriminate electrons from background. The 
background, which is mainly due to pions misidentified as electrons and to a lesser extent due 
to kaons and anti-protons, is largely suppressed by combining the two independent estimators 
into a combined estimator D e \ e , as explained in appendix [A] The three estimators are defined 
such that D = 1 for genuine electrons and D = for pion background. 

The calorimeter based electron identification [47] is track seeded, which means the cluster 
shape estimators are calculated from energy deposits in LAr calorimeter cells lying within a 
cylinder of 30 cm around the extrapolated track trajectory. The Cluster energies are corrected 
for energy losses in the dead material in front of the LAr calorimeter. Electron candidates with 
energy depositions close to inactive regions between LAr calorimeter modules are rejected. 
Five estimators are defined: four cluster shape variables and the ratio of the energy deposited in 
the electromagnetic part of the LAr calorimeter to the momentum of the corresponding track. 
These estimators together with the logarithm of the total energy and the z position of the clus- 
ter, are mapped onto one single estimator -D ca i using the artificial neural network Multilayer 
Perceptron [|50l . 

The measured specific ionisation loss of the track, dE/dx, is translated into x 2 -probabilities 
of corresponding particle hypotheses P (dE/dx, e) for electrons and P '(dE/dx, it) for pions, 
which constitute the main background. From both probabilities the estimator 

P(dE/dx,e) 

dE/dx p( dE / dXj e ) + p^E/dx, 7T) ^ ' 

is constructed. The simulation of the specific ionisation was studied in detail in order to describe 
precisely the measured energy losses fl5Tl . 

The performance of both discriminator variables is validated using Monte Carlo and data 
samples of identified electrons and pions in the transverse momentum range of interest, 1 < 
Pr(e) < 5 GeV, selected in decays J/ij) — > e + e~ and K® — > by means of the "tag and 

probe method" B71 . 

The simulation describes well the distribution of the discriminators -D ca io and D d E/dx as 
measured in data, as can be seen in figure |2] The deviations of the simulation from the data 
at small D values in the electron sample are due to a small remaining pion contamination in 
the data. Also the combined estimator D ele is found to be well described by the simulation and 
shows an excellent separation of the electron signal from the pion background. Isolated elec- 
trons are selected for D ele = 0.825 with an efficiency of more than 90% for a pion background 
rejection of about 99%. 

4.3 Event Selection 

A di-electron sample is obtained by selecting events with two or more offline reconstructed 
electron candidates, requiring D e \ e > 0.825. To account for the P T resolution of the third 
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trigger level, the cut on electron tracks reconstructed offline is raised by 100 MeV above 
the trigger threshold of the respective subtrigger (see tabled)), which recorded the event. This 
implies two electrons with P T > 1.3 GeV for the low-Py subtrigger and one electron with 
P T > 1.6 GeV (P T > 2.1 GeV) for the medium- P T (high-P T ) subtrigger. 

In order to remove background from non-ep sources, the measurement of a good event 
vertex is required. The event vertex is reconstructed from all charged tracks of an event and 
its position along the beamline has to be within 30 cm around the nominal interaction point. In 
addition, timing vetoes are applied to further reduce non-ep interaction induced backgrounds. 

In order to reject background from DIS, events with a positron in the LAr calorimeter iden- 
tified by the standard electron identification [1531 and with E(e + ) > 8 GeV are rejected. As 
the Pt{&) -distribution of semileptonically decaying 6-quarks falls steeply, almost all 6-decay 
positrons are at low energies and thus not affected by this cut. Also DIS events are rejected 
which have an electromagnetic cluster in the SpaCal calorimeter with energy above 8 GeV con- 
sistent with originating from the scattered beam positron. Events with Q 2 < 2 GeV 2 are not 
rejected by these cuts, since the beam positrons leave the detector undetected along the beam 
pipe. 

Only events with measured inelasticities in the phase space region of this measurement, 
0.05 < yh < 0.65, are accepted. The inelasticity variable is reconstructed from the sum over 
all final state particles y h = Yl/ii^i ~ Pz,i) / (^E e + heaxa ), where -Ee+beam denotes the energy 
of the beam positron. Particles belonging to the hadronic final state (HFS) are reconstructed 
using a combination of tracks and calorimeter deposits in an energy flow algorithm that avoids 
double counting [54]. Ei and P Zji denote the energies and longitudinal momenta of all final state 
particles, which correspond to the visible hadronic final state in case of photoproduction, and in 
case of DIS background also includes the scattered positron. The upper cut on the inelasticity 
suppresses effectively remaining DIS events. 

The beauty signal is further enriched by rejecting electron candidates, which are in a dense 
hadronic environment. For this purpose the variable -RE.cone is defined as the ratio of the summed 
energy of all HFS particles in a cone of 18° around the electron track direction, E conc , to the 
electron energy E e , which must not exceed an upper threshold: 

P E ,cone = ^ < 350% . (3) 

The effect of this cut is twofold: First, it reduces misidentified electron candidates resulting from 
overlapping showers in the LAr calorimeter. Second, it enriches electrons from semileptonic 
beauty decays, which are in general isolated from hadrons due to the large 6-mass. 

Finally, electrons from photon conversions are rejected by the three following cuts. First, 
the distance of closest approach in the transverse plane, dca e , of the electron tracks to the beam 
line is restricted to be smaller than 0.2 cm. Second, a photon conversion finder searching for 
displaced vertices is used to identify electrons originating from the photon conversion process 
7 — > e + e~. Third, the invariant mass of the selected electron pairs is required to be m e i e2 > 
1.2 GeV. This cut rejects e + e _ pairs from Dalitz decays and most of the remaining background 
from photon conversions. 

The selection cuts are summarised in table |2] After applying all cuts about 1500 electron 
pairs are selected. In the rare cases with more than two selected electrons per event all pair 
combinations are considered in the analysis. 
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Overview of the Selection Cuts 

Trigger selection 

• track multiplicity cuts 

• 1 or 2 online identified electrons 
Offline electron selection 

• 2 electron candidates with: 

- D ele > 0.825, R E>cone < 350% 

- P T (e) > 1 GeV, 20° < 0(e) < 140° 

- verification of the L3 P T (e) -thresholds 100 MeV above the P T (e) -threshold of the 
respective subtrigger which recorded the event (see tabled] and text) 

Background rejection and further cuts 
Rejection of non ep-background: 

• good vertex, timing vetoes 
Rejection of DIS events: 

• no identified scattered beam positron 

• 0.05 <y h < 0.65 

Rejection of photon conversions and Dalitz decays: 

• rn e i t e2 > 1.2 GeV 

• no converted photon 

• dca e < 0.2 cm 

Table 2: Overview of the online and offline selection cuts. More details on the selection proce- 
dure can be found in Wf\ . 

5 Data Analysis 

The selected di-electron sample is dominated by events from inelastic J/^-meson production. 
While decays of J/^-mesons can be easily identified by kinematic reconstruction of the J/ij; 
mass peak, the separation of the bb signal events from the other backgrounds originating from 
the production of light quarks and open charm production is more difficult. 

In the following, the reconstruction of the transverse momentum of the produced 6-quarks 
and the flavour separation of the different processes are described. 

5.1 Reconstruction of 6-quarks 

The transverse momentum of 6-quarks is reconstructed for the measurement of the differential 
cross section da/d(PT(b)), where (P T (b)) is the quadratically averaged transverse momentum 
of the b and b quark as defined in equation [T] The bb cross section is largest at small transverse 
momentum at (P T (b)) « m&, a kinematic region where standard jet finders cannot be used due 
to isotropic decay topologies. Therefore an alternative, referenced as the hemisphere method, 
is exploited. This method was applied in a previous analysis [[56*1 to reconstruct the directions 
and momenta of charm quarks in the production of cc-pairs in DIS, and is also well suited to 
reconstruct the transverse momenta of 6-quarks in bb production ll47l . 



11 



As illustrated in figure |3] an event is divided into hemispheres, using the thrust-axis which 
is calculated in the laboratory frame in the plane transverse to the beam directions (x-y plane). 
Using the transverse momenta from all particles of the HFS, the thrust-axis in the transverse 
plane is given by the vector a maximising the sum of the projected transverse momenta onto it, 

T = m a x(a)( ^ UFslS 'f T A with|a| = l . (4) 

\ SieHFS \ P T,i\ J 

A plane perpendicular to the thrust-axis defines two hemispheres, one of them containing the 
fragmentation products of the 6-quark, and the other one containing the fragmentation products 
of the 6-quark. 

Two observables -P^hem.i an d -Pr,hem.n are used to reconstruct the mean transverse momen- 
tum of the b (b) quark produced in the hard interaction. These observables, which are derived 
from the HFS particles assigned to the corresponding hemispheres, show a good correlation 
to the transverse momentum of the b (b) quarks in the hard process. However, the hadronic 
final state also contains particles from the so called proton remnant, leaving the interaction in 
the positive ^-direction of the detector and thus deteriorating the above correlation. Simulation 
studies show that the correlation with the 6-quark transverse momentum is improved by exclud- 
ing particles in the forward direction at polar angles below 15 degrees. The transverse momenta 
of the b (fr)-quarks are therefore approximated by: 

^T.hem.I (hem.IT) = ^T,i with ^ > 15° . (5) 

i £ hem. I 
(i £ hem. II) 

This reconstruction method is very reliable at large Px,b, where two hard jets are measured 
in the final state. At small Pr(b) the transverse momenta of HFS particles in the hemispheres 
are mainly generated by the b and 6-hadron decays themselves and are related to the mass of 
the 6-quark: l-Pr^em.il ~ |-Pr,hem.n| ~ m b- In order to allow for a good reconstruction of Pr(b) 
down to the bb production threshold, i.e. Pr(b) ~ GeV, the average transverse beauty mass is 
used: 

mr(b) = yjml + (P T (6)) 2 . (6) 

Detailed studies [47] demonstrated that the average transverse beauty mass can be well recon- 
structed from the experimental observables -Pr.hem.i and -Pr,hem.n using the relation: 

m T ,rec(b) = a* (|Pr,hem.l| + l-P^hem-Ill)/^ , (7) 

with a being a constant parameter set to a = 1.09, such that the correlation between generated 
and reconstructed m T {b) is maximised. This correlation as obtained by simulation is shown in 
figure SI For values of m& it in the range 4.5 < m b < 5.0 GeV the dependence of this correlation 
on m b is negligible. 
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j->max(el,e2) 
-^EjCone 


"^ele 




0.825 - 0.875 


0.875 - 1.0 


150 - 350% 


Bl 




- 150% 


B2 


S 



Table 3: Definition of the four regions used to constrain mainly the uds background. 



5.2 Quark Flavour Separation 

For the discrimination of the bb signal against remaining background from misidentified elec- 
trons and for the separation of the different quark-flavour components contributing to the di- 
electron signature, a template method is used. Several independent phase space regions are 
defined such that individual background sources are enhanced in certain regions of the phase 
space and can be tested while other contributions are suppressed. Finally the 6-signal ("beauty") 
and the background contributions are obtained by an unfolding procedure. Background sources 
determined by this method are the production of light quarks ("uds"), open charm production 
("charm") and the production of J/^-mesons ("J/^>"). The uds background contains also a 
small fraction of charm and beauty events, where at least one electron candidate does not origi- 
nate from a semi-electronic heavy quark decay. 

5.2.1 Fraction of Light Quarks 

In order to determine the background contributions due to misidentified electrons the data sam- 
ple is grouped in four regions Bl, B2, B3 and S using different electron quality criteria on 

r - 1 min(el,e2) , 7 -,max(el,e2) , , , n,! ,-.min(el,e2) < T-,max(el,e2) • • , 

D ele and R Eco ly e , see table |3J D cle and -R Ec one 316 tne minimum and maxi- 
mum value of Dele and -fircone respectively, of the two electron candidates, which form the elec- 
tron pair. Bl, B2 and B3 are background enhanced regions and S denotes the electron signal 
enhanced region, which is defined by tight electron identification and isolation cuts. Templates 
for the determination of the background fractions are obtained from Monte Carlo simulations, 
see figure IU More than 70% of the beauty, charm and J/ip events populate the signal enhanced 
bin S, since these events contain genuine electrons. The uds events are enriched in the three 
background bins Bl, B2 and B3, due to misidentified electrons. The measured number of 
events in these three background bins constrain mainly the uds background fraction. 

5.2.2 Heavy Quark Fractions 

In the signal enhanced region S, the individual contributions from beauty, charm and J/ip can be 
disentangled by investigating the charge product, q el ■ q e2 , of the e ± -candidates, their azimuthal 
separation A0 el e2 = |0 el — (p e2 \, and their invariant mass m el e2 . Templates of the different 
background sources and of the beauty signal, which are all restricted to the signal enhanced 
region S, are shown in figure |6]as function of the invariant mass m e i,e2 and the signed azimuthal 
separation A0 e i je2 • q el ■ q e2 . 
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Phase Space 

Q 2 < 1 GeV 2 
0.05 < y < 0.65 
M6)|, M6)|< 2 

Table 4: Definition of the kinematic range of this measurement. 



The different templates show specific characteristics: J /ip events have oppositely charged 
electrons and cluster at m e i e2 = mj/^, whereas background from open charm production cov- 
ers a large mass range. Electrons from open charm decays are found mostly back-to-back and 
with opposite charge sign, whereas electron pairs from beauty decays populate all A0 el e2 val- 
ues with both charge sign combinations. Both charge products are also found in the uds back- 
ground, which however populate on average regions with smaller m e i e2 values. Large values 
of m e i te 2 are solely populated by beauty decays. 

These distinct signatures of the individual background sources, i.e. uds, J ftp — >• ee and 
cc ee, are exploited by dividing the signal enhanced region 5 into 12 subregions (51 to 512) 
as shown in figure |6] In the following the three background enhanced bins B1-B3 and the 12 
signal enhanced bins 51-512 are referred to as "Flavour Separator", for which templates are 
derived. 

5.3 Unfolding 

Using an unfolding procedure the number of background events iV u ds, Nj/^, A^harm and the 
number of beauty events Abeauty,i in four bins of (Pt(&)) are derived. A regularized unfolding 
procedure is used with a smoothness condition. The procedure is explained in appendix [Bj All 
efficiency corrections and migration effects are described by the response matrix A, which cor- 
relates the number of reconstructed events in the Flavour Separator distribution in bins of m T rec , 
represented by the vector y, with the distribution x on parton level via the matrix equation 

y = Ax + b . (8) 

The vector x, defined as x T = (xb eauty , £ c harm, X J/^-, ^uds), contains contributions from beauty 
binned in (Pr(&)), charm, J/tp and uds. The contribution from beauty (xb eau ty) is defined 
according to the phase space given in table |4] The vector b contains the background contribu- 
tion from DIS events, which is taken from simulation. All other background contributions are 
incorporated in the response matrix and are determined by unfolding. 

Signal and background templates as function of nine m T rcc bins are generated by Monte 
Carlo simulations and fitted to the data. The unfolding procedure uses in total N mTicc x 
^Flavour Separator = 9 x 15 input bins and determines the three background fractions and the 
number of beauty events iV beautyj j m four (P T (b)) bins. A schematic representation of the proce- 
dure is shown in figure [7] In this procedure the m T rcc dependences of the different background 
contributions from uds, J/ip — > ee and cc — > ee are fixed by the Monte Carlo predictions. The 
latter is motivated by recent measurements of the differential cross sections of charm production 
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Background Correlations 


Fractions 




beauty charm J/ip uds 


[%] 


beauty 


1 -0.46 -0.18 -0.18 


25.8 ±3.6 


charm 


1 -0.03 -0.27 


17.6 ±3.3 


J/i> 


1 0.03 


29.0 ±2.1 


uds 


1 


25.3 ±3.0 



Table 5: Correlations between the signal (beauty) and the different background contributions, 
and the determined relative fractions with their errors for the data sample. The fraction of DIS 
events (not given in the table) is 2.3%. 



at HERA, which were found to be consistent with theoretical models and Monte Carlo programs 
used in this analysis H57U58H . 

The fitted beauty signal and background contributions are shown in figure [8] in the three 
background and in the signal enhanced regions. The event numbers resulting from the fitted 
fractions show very good agreement with the data considering statistical errors only. A clear 
enhancement of the genuine electron signal due to the tightening of the electron identification 
cuts is seen when going from the first background enhanced region (BY), which contains more 
than 80% uds background to the signal enhanced region (51-512) with less than 20% of uds 
background. 

The correlations between the beauty signal and the background sources, which are largest 
between beauty and charm, are given in table |5] together with the determined fractions of the 
selected data sample. 

The distribution of the data as a function of the Flavour Separator is shown in figure [9] 
together with the result from the fit of the beauty and the various background contributions. 
Good agreement is found considering statistical errors only. 

Control distributions of electron variables are presented in figure QJJfor the electron enriched 
signal region (51-512). The data are compared to the simulated beauty signal and background 
distributions using the quark flavour decomposition determined by the unfolding procedure. 
The main characteristics of the signed variables A</> el e2 ■ q e \ ■ q e 2 and m e i je2 ■ g el ■ q e2 , and 
-Pr(e) and $(e) are well described by the Monte Carlo simulation. In figure [TT| additional 
control distributions are presented showing the P^-spectra of the three highest P^-tracks. These 
distributions are strongly dependent on the track trigger conditions used, and imperfections of 
the trigger simulation would be visible here. 

Reasonable agreement between the data and the Monte Carlo simulation is obtained in all 
distributions which gives confidence that the Monte Carlo simulation is able to correctly model 
the detector response used for the unfolding procedure. 
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5.4 Cross Section Determination and Systematic Uncertainties 

The visible cross section is measured for the phase space as defined in tablelU The bin-averaged 
differential cross section is obtained as 

da(ep -> ebbX) _ N heauty ,i 

d(P T (b)) ~ C-BR-A{P T>i (b)) ' W 

where C is the luminosity, A(Pr,i(b)) the bin width, iVb eau ty,i the number of unfolded signal 
events in the corresponding bin and BR = 6.17% the effective branching fraction computed 
from [55] for a bb pair decaying into at least two electrons. For the calculation of cross section 
uncertainties correlations between bins are taken into account. 

The systematic uncertainties related to the measurement of the number of bb signal events 
are listed in the following. The effect on iV beauty ; is calculated by varying the sources of uncer- 
tainties in the simulation and by propagating these variations to the measurement through the 
response matrix A and the background term b in equation [8] 

• The uncertainty on the electron identification is determined using J/ijj — > e + e~ events 
(see figure |2]), by comparing the distributions of the electron discriminator Z) ele between 
data and Monte Carlo around each of the used cut values of D c \ c = 0.875 and D c \ c = 
0.825. The cut on D e \ e is varied in MC by ±0.025 which covers any possible shift in the 
-D e i e distribution between data and simulation. This cut variation on D e \ e propagated to 
the total beauty cross section results in an uncertainty of ±6.8%. 

• The uncertainty on the track finding efficiency of electrons is conservatively estimated to 
be 2% per track resulting in an uncertainty of the total beauty cross section of ±4%. 

• The trigger uncertainty of the FTT at levels LI and L2 are about 1 — 2% each. The dom- 
inating contribution to the trigger uncertainty is due to the uncertainty of the calibration 
constants of the JT used at L3. To quantify this uncertainty, the JT calibration constants 
used in the simulation are varied by scaling the default calibration constants by 15% B7l . 
The systematic error on the total beauty cross section due to the uncertainty on the trigger 
efficiency is determined to be ±8.6%. 

• Model uncertainties of the beauty signal are determined by comparing the default re- 
sponse matrix computed by taking the average of the two Monte Carlo samples (CAS- 
CADE and PYTHIA) with two alternative response matrices based on one of the Monte 
Carlo samples. The relative maximum difference with respect to the default response ma- 
trix is computed for each entry of the matrix and propagated to a model uncertainty on 
the total beauty cross section of ±3.3%. 

• The uncertainty of the charm contribution is evaluated from the relative difference be- 
tween the Monte Carlo generators CASCADE and PYTHIA in a similar way as for the 
beauty signal. The systematic error on the extracted total beauty cross section due to the 
charm model is determined to be ±3.6%. 



16 



• The uncertainty due to the fragmentation function of the heavy quarks is estimated by 
reweighting the events according to the longitudinal string momentum fraction z carried 
by the heavy hadron in the Lund model using weights of (1 =F 0.7) • (1 — z) + z • (1 ± 0.7) 
for charm quarks and by (1 =p 0.5) • (1 — z) + z - (1 ± 0.5) for beauty quarks [|2~2l . 
The corresponding systematic error on the total beauty cross section is determined to be 
±3.4% resulting from charm and ±2.2% from beauty fragmentation uncertainty. 

• The uncertainty on the contribution from the remaining uds background due to misidenti- 
fied and real electrons, was determined by varying their relative contributions by a factor 
two up and down. The corresponding systematic error on the total beauty cross section is 
determined to be ±3.4%. 

• CASCADE does not fully simulate the radiative tail of J ftp — > ee events. To estimate the 
uncertainty on the modelling of it, weights are applied, which are obtained from an elastic 
J/ip — > ee simulation with radiative QED corrections [|59l . The systematic uncertainty 
on the total beauty cross section is estimated to be ±3.5%. 

• The uncertainty of the DIS-background, represented in equation [8] by the vector b, is 
taken to be 100% and results in an error on the total beauty cross section of ±4.5%. 

In addition, a global normalisation uncertainty of 4.1% is included with contributions from 
the integrated luminosity uncertainty of 2.7% and from the uncertainty on the semi-electronic 
branching fractions of 3.0%. 

Adding all above contributions in quadrature gives a total systematic error of 15.4% on the 
total beauty cross section. 



6 Results 

The differential cross section da(ep — > e bbX) /d(Pr(6)) is measured in the phase space defined 
in table H] using the unfolding procedure as described in section [531 

The result is shown in table [6] together with statistical and total errors and the coefficients 
describing the statistical correlations between bins. In order to cross check the unfolding pro- 
cedure the cross section extraction is repeated without regularisation condition. The results 
obtained with and without regularisation are found to be consistent within the uncertainties. 

The measured differential beauty cross section is compared in figure [12] with an NLO QCD 
prediction in the fixed flavour number scheme as calculated by the program FMNR. The figure 
also shows the ratio of the measured cross section and the NLO QCD cross section. The uncer- 
tainties of the measurement are smallest at low (Prib)), where the cross section is largest. The 
theoretical prediction of the differential cross section agrees with the measurement within the 
large experimental and theoretical uncertainties. The prediction has a tendency to be below the 
data, a trend also observed in previous beauty cross section measurements at large transverse 
beauty momenta. 
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By integrating the differential cross section the total inclusive beauty photoproduction cross 
section is measured as: 

a(ep ->■ e bbX) = 3.79 ± 0.53 (stat.) ± 0.58 (sys.) nb , (10) 

to be compared with the NLO prediction obtained from FMNR of <r(ep — > e bbX) = 2.40^^1 nb. 
The measured cross section is higher, but within the large experimental and theoretical uncer- 
tainty consistent with the NLO expectation. 



7 Conclusions 

The inclusive and differential cross section of beauty photoproduction was measured in the di- 
electron final state, using the HI detector at the HERA collider. The cross section is measured 
as function of the quadratically averaged transverse momentum of the produced beauty quarks 
(Pr(b)), with a special focus on the low (P T (b)) regime. Background from uds, charm and J/ip 
production is determined exploiting angular, charge and mass correlations of electron pairs in 
an unfolding procedure. 

The measured cross section is compared to a QCD prediction at NLO performed in the fixed 
flavour number scheme and evaluated with /j, R — /j, F — 1/2 \/m\ + (P T {b)) 2 as choice for the 
renormalisation and factorisation scale. The NLO prediction lies below the data but within the 
large experimental and theoretical uncertainty they agree. 

This measurement is in good agreement with previous measurements of beauty photopro- 
duction at HERA and it extends the previously experimentally accessible phase space towards 
the beauty production threshold. 
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A Electron Discriminator Combination 



The track seeded and calorimeter based electron discriminator -D ca io and the electron discrim- 
inator DdE/dx, based on the specific energy loss measurement in the CTD, are mapped to a 
combined discriminator using the expression 



D e ie (a, b, c, d) 



'l-D calo ) a -l\ c -\{l-D dE/dx ) b -l\ 



((1 - D^Y - 1) ■ ((1 - D dE/dx ) b - 1) + (1 - D calo )- ■ (1 - D dE/dx f ' 

(11) 



which for the parameter choice a — b — c — d — 1 corresponds to Bayes' theorem: 

D ele (l, 1, 1, 1) = 3 \D^\D dE/dx \ (12) 

-t^calo • -^dE/dz + U — -t^caloj • [i- ~ ^dE/dx) 

However, in order to obtain a sensible mapping behaviour of -D calo and D dE i dx onto D c \ c when 
their respective values are close to 1 and or both of them are close to 1, the parameters 
a = b = 0.6 and c = d = 1.05 are chosen. 



B Unfolding procedure 

The differential cross section of beauty photoproduction is extracted from the measured di- 
electron spectrum using an unfolding procedure as implemented in TUnfold fl6TTl . 

The vector y, representing the number of measured events, is related via the matrix equation 
y = A ■ x + b to the true distribution represented by a vector x, which is determined by 
unfolding. The response matrix A describes the detector acceptance, contains all selection 
efficiencies and takes migration effects between bins into account. Additional background, not 
determined by the unfolding procedure and taken from external information, is represented by 
the vector b. 

An estimator x of the true distribution x is obtained by unfolding the measured distribu- 
tion y. For the construction of x additional assumptions, e.g. on the smoothness of the de- 
convoluted distribution (regularisation), and an additional constraint on the number of observed 
events are applied. In general x is obtained by minimising a x 2 function given by: 

X 2 (x,r,/i) :=Xa(x) + t-x£(x)+//-Xat(x) . (13) 

This equation describes the minimisation of the unfolding problem xi( x ) with the two side 
conditions given by an d Xjv(x)- 

The actual minimisation problem is defined by the standard x 2 function: 

X 2 A m := \ (y - b - Axf V- 1 (y - b - Ax) (14) 

with V = cov(yi,yj) being the covariance matrix of the data. This function minimises the 
deviation of the estimator Ax from the measured, and background subtracted vector y — b. 
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The additional constraints are given by: 



Xi(x) := x T Lx 
X%{x) := n ohs - 



) 



(15) 



(16) 



with L being the regularisation matrix, m the number of reconstructed bins and n ^ s the total 
number of observed events after background subtraction, which ensures that the total number 
of events is conserved. Both functions enter equation [T3l with the parameters r and /i, where r 
is often denoted as regularisation parameter and fi as Lagrange Multiplier. 

The x!( x ) function is a measure for the smoothness of the result. The matrix L is chosen 
such that the second derivative of x between bins describing beauty production is minimised. 
The regularisation parameter r determines the strength of the smoothness constraint. For the 
regularised unfolding r is chosen such that the correlations of the covariance matrix of the 
unfolded distribution x are minimised 11621 . 
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HI Beauty Photoproduction Cross Sections 
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Table 6: Differential cross sections for the phase space defined in table |4] obtained from unfolding with regularisation condition. Also given 
are the statistical and total errors, the coefficients of the statistical correlations, the used bin boundaries in (Pr(6)) and the corresponding bin 
centres [|60l {Pr,bc{b)). The remaining columns list the the bin-to-bin correlated systematic uncertainties in the cross section measurement 
due to uncertainties of the beauty (8 b s ) and charm (<^ s ) modelling, the uds background (5^), the electron identification the 

beauty (5^ s b ) and charm (5/^ c ) fragmentation, the modelling of the radiative tail of J/tp — > e + e~ events (Ssyf.) and the DIS background 
((*)£i S )- Not listed in the table is the 4.1% normalisation uncertainty. 
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Figure 2: Normalized discriminator distributions for the separation of electrons and pions as 
obtained from J/ip — > e + e~ and — > decays using the tag and probe method, a) the 

track seeded, calorimeter based discriminator -D ca i , b) the discriminator D dE / dx based on the 
measurement of the specific energy loss in the CTD and c) their combination D e \ e . Data are 
represented by circles and Monte Carlo simulations by histograms. 
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A 



thrust 

x 

Figure 3: Schematic illustration of the determination of the thrust axis in the plane transverse 
to the ep beams. The transverse thrust axis, indicated by the dashed arrow, maximizes the sum 
of momenta projected onto it in this plane. The thrust axis allows the event to be divided into 
two hemispheres, each containing the decay products of a beauty quark, used to reconstruct the 
average transverse beauty mass m T ,iec(b) as defined in equation [7J 




<m T (b)> [GeV] 

Figure 4: Correlation between the reconstructed transverse beauty mass mr,rec(&) and the trans- 
verse mass {rriT (&)) calculated from the quadratically averaged transverse momentum of the 
generated beauty quarks. The inner line on the diagonal indicates the correlation of mT, T ec(b) 
and (rriT(b)), and the outer two lines show the la error band. The used binning (dotted grey 
lines) for the vectors x and y entering the unfolding procedure are also shown. 
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Figure 5: Templates used to separate the light quarks (uds) from the heavy quark flavours as 
obtained by the Monte Carlo simulation. For the definition of the background enhanced regions 
B1-B3 and the signal enhanced region S see tableland text. 
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Figure 6: Templates in the plane A0 elie2 • q e \ ■ q e2 , m e \ fi 2 and restricted to the signal enhanced 
region S used to separate the heavy quark flavours as obtained by the Monte Carlo simulation. 
Also shown is the bin numbering S1-S12 of the 12 subregions of S. For the definition of the 
signal enhanced region S see tableland text. The two vertical lines indicate the peak invariant 
mass region of the J/ip — >■ e + e~ decays. 
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Figure 7: Structure of the response matrix A used to solve the matrix equation y = A • x + b by 
unfolding. The indicated numbers specify the number of used bins. The sub-binning in y given 
by the Flavour Separator allows the discrimination of the beauty signal from the uds, charm and 
J ftp backgrounds. See text for details. 
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Figure 8: Number of di-electron events in the background and signal enhanced regions as de- 
fined in table S Data are represented as points with the statistical error indicated by the error 
bars. Also shown in colour is the decompositions of the event yields as determined by the 
unfolding procedure. 
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Figure 9: Number of di-electron events in the flavour separator histogram compared to the num- 
ber of fitted events and their decomposition. Data are represented as points with the statistical 
uncertainties indicated by the error bars. The bin numbering scheme as defined in figure [6] and 
table [3] is used. 
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Figure 10: Control distributions of the electron candidates compared to Monte Carlo simula- 
tions using the quark flavour decomposition determined by the unfolding procedure: a) signed 
azimuthal separation A0 el e2 ■ q e \ ■ q e 2 defined by the charges multiplied with the azimuthal angle 
difference of the two electron candidates, b) signed invariant mass m e i e2 • q e \ • q e 2 defined by 
the charges multiplied with the invariant mass of the two electron candidates, c) polar angle of 
the electron candidates and d) transverse momentum of the electron candidates. Data are repre- 
sented as points with the statistical uncertainties indicated by the error bars. The distributions 
are restricted to the electron enriched region (S). 
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Figure 1 1 : Control distributions for the three highest P T -tracks of the hadronic final state as 
function of the track P T . Data are compared to the Monte Carlo simulations using the the quark 
flavour decomposition determined by the unfolding procedure. Data are represented as points 
with the statistical uncertainties indicated by the error bars. 
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Figure 12: Differential beauty cross section da / cI(Pt(&)) shown as function of the quadrati- 
cally averaged transverse momentum of the beauty quarks (P T (b)) (upper part). The data are 
represented by points with inner vertical error bars representing the statistical errors and outer 
error bars representing the total error. The vertical gray lines indicate the bin boundaries in 
(P T (b)) of each data point and the points are shown at the bin centred positions. The data are 
compared to the FMNR NLO QCD calculation (solid line) with the uncertainty represented as 
shaded band. Also shown is the ratio of the measured cross section to the calculated NLO QCD 
prediction, ^—f (lowerpart). 
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